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NiCl2-4SC(NH2)2 (DTN) is a quantum S = 1 chain system with strong easy-pane anisotropy and a new 
candidate for the Bose-Einstein condensation of the spin degrees of freedom. ESR studies of magnetic exci- 
tations in DTN in fields up to 25 T are presented. Based on analysis of the single-magnon excitation mode in 
the high-field spin-polarized phase and previous experimental results [Phys. Rev. Lett. 96, 077204 (2006)], 
a revised set of spin-Hamiltonian parameters is obtained. Our results yield D = 8.9 K, J c — 2.2 K, and 
Ja.b = 0.18 K for the anisotropy, intrachain, and interchain exchange interactions, respectively. These values 
are used to calculate the antiferromagnetic phase boundary, magnetization and the frequency-field dependence 
of two-magnon bound-state excitations predicted by theory and observed in DTN for the first time. Excellent 
quantitative agreement with experimental data is obtained. 
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Antiferromagnetic (AFM) quantum spin-1 chains have 
been the subject of intensive theoretical and experimental 
studies, fostered especially by the Haldane conjecture []]]. 
Due to quantum fluctuations, an isotropic spin-1 chain has a 
spin-singlet ground state separated from the first excited state 
by a gap A ~ 0.41J J2], where J is the exchange interac- 
tion. As shown by Golinelli et al. J^, the presence of a strong 
easy-plane anisotropy D can significantly modify the excita- 
tion spectrum, so that the gap size is not determined by the 
strength of the AFM quantum fluctuations exclusively, but 
depends on the dimensionless parameter p = D/J. The Hal- 
dane phase is predicted to survive up to p c = 0.93 J3l, where 
the system undergoes a quantum phase transition. For p > p c 
the gap reopens, but its origin is dominated by the anisotropy 
D, and the system is in the so-called large-D regime. While 
the underlying physics of Haldane chains is fairly well under- 
stood, relatively little is known about the magnetic properties 
(and particularly the elementary excitation spectrum) of non- 
Haldane 5=1 AFM chains in the large-D phase. Intense the- 
oretical work and numerous predictions il BBS Hi 
make the experimental investigation of large-Z) spin-1 chains 
a topical problem in low-dimensional magnetism. 

Recently, weakly-coupled spin-1 chains have attracted re- 
newed interest due to their possible relevance to the field- 
induced Bose-Einstein condensation (BEC) of magnons. 
When the field H, applied perpendicular to the easy plane, 
exceeds a critical value H c \ (defined at T = 0), the gap closes 
and the system undergoes a transition into an XF-like AFM 
phase with a finite magnetization and AFM magnon excita- 
tions. If the spin Hamiltonian has axial symmetry with re- 
spect to the applied field, the AFM ordering can be described 
as BEC of magnons by mapping the spin-1 system into a gas 
of semi-hard-core bosons itllll . The applied field plays the role 



of a chemical potential, changing the boson population. In ac- 
cordance with mean- field BEC theory ||12lll3Lll4ll, the phase- 
diagram boundary for a three-dimensional system should obey 

3/2 

a power-law dependence, H — H c i ~ T c ' , as T — > 0. Above 
the upper critical field, H C 2 (defined at T = 0), the system is 
in a fully spin-polarized (FSP) phase, and the excitation spec- 
trum is formed by gapped magnons. 

The compound NiCl2-4SC(NH 2 )2 (dichloro-tetrakis 
thiourea-nickel(II), known as DTN), is a new candidate 
for studying the BEC phenomenon in magnetic fields. It 
has a tetragonal crystal structure, space group 14, with two 
molecules in the unit cell 1 15, 3. The anisotropy, intra- and 
inter-chain exchange parameters, D = 8.12 K, J c = 1.74 K, 
and J a .b = 0.17 K, respectively, were obtained from a fit of 
zero-field inelastic-neutron-scattering data (l7\\ using gener- 
alized spin- wave theory 1 1811 . It was shown that the Ni spins 
are strongly coupled along the tetragonal axis, making DTN 
a system of weakly interacting 5=1 chains with single-ion 
anisotropy larger than the intra-chain exchange coupling. 
Recently, it was proposed flril 17 1 that the field-induced 
low-temperature transition of DTN to the AFM-ordered state 
can be interpreted as a BEC of magnons, with values of 
H c i = 2.1 T and H C 2 = 12.6 T for the lower and upper 
critical fields, respectively. Although the observed overall 
picture is consistent with the BEC scenario, some questions 
remain. One of the unsolved problems is a pronounced 
(about 2 T) disagreement between the value for the upper 
critical field obtained experimentally and the predicted one 
(H C 2 = 10.85 T) U3]. This discrepancy is rather puzzling, 
particularly in light of the excellent agreement for H c \. 

In this Letter, we report electron spin resonance (ESR) stud- 
ies of the elementary excitation spectrum in DTN in magnetic 
fields up to 25 T (~ 2H C 2). A distinct advantage of the high- 
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FIG. 1: (color online) A schematic view of the magnetic-excitation 
dispersion in an 5 = 1 Heisenberg chain with strong easy-plane 
{D > 0) anisotropy for two arbitrary fields, H < H c i (left) and 
H > H C 2 (right), with a magnetic field applied along the principal 
axis z. Note that the ESR transitions denoted by A, B, C, E, and F 
occur at k = 0. Two-particle continua predicted for both regions are 
not shown for simplicity. 



field approach for determining the spin-Hamiltonian parame- 
ters of spin-1 AFM chains is the availability of exact theo- 
retical expressions for the spin-polarized phase. Analysis of 
the high-field single-magnon branch has allowed us to obtain 
the precise value for the anisotropy parameter, D = 8.9 K. 
Using results of recent thermodynamic and neutron-scattering 
measurements JlH. 17], we were able to refine the exchange 
parameters, obtaining J c = 2.2 K, J a ^ = 0.18 K for intra- 
and inter-chain exchange interactions, respectively. In addi- 
tion, we present magnetocaloric-effect and low-temperature 
magnetization data allowing us to check the obtained set of 
parameters. These values agree well with those obtained by 
fitting the AFM phase boundary and the magnetization with 
results of quantum Monte Carlo simulations, and nicely re- 
produce both critical fields, H c \ = 2.1 T and H C 2 = 12.6 T. 
Furthermore, we report on the first direct and reliable obser- 
vation of the two-magnon bound states in spin- 1 AFM chain 
system with strong easy-plane anisotropy, predicted by theory 
J9D for the high-field spin-polarized phase. Using the obtained 
parameters of the spin-Hamiltonian we were able to calculate 
frequency-field dependences of the two-magnon bound state 
excitations. Excellent quantitative agreement between the the- 
oretical predictions and experiment is obtained. 



The investigation of the magnetic excitation spectra in 
single-crystalline DTN samples was done using a tunable- 
frequency submillimeter-wave ESR spectrometer |[l9tl with 
external field H applied along the tetragonal c axis. 

In sufficiently small fields, a spin-1 Ni(II) system is in the 
quantum-paramagnetic (QPM) phase, having the M s = 
ground state. The excitation spectrum is gapped and deter- 
mined by the AM S = ±1 transitions. For fields H applied 
perpendicular to the easy plane, the Hamiltonian can be writ- 



n = J - S i ■ s j+^ + J2i D ( s D 2 + mbHS?], (d 

where v = {a, b, c}. A corresponding dispersion of low-field 
magnetic excitations is shown schematically in Fig. Q] (left). 
Then, the frequency-field dependence of the AM S = ±1 
ESR transitions can be written as uj a / b = A ± g/iBH, 
where A corresponds to the energy gap at k = and H = 0. 
These transitions were observed in our experiments and are 
denoted in Fig.|2]by circles I20I1 . The zero-field splitting mea- 
sured directly yields A = 269 GHz, which agrees quite well 
with results of neutron-scattering measurements, Ak=o = 1.1 
meV fnil . No splitting within the excited doublet in zero 
field was observed, which indicates the absence of in-plane 
anisotropy. From the slope of the frequency-field dependences 
of the branches A and B, the g factor was determined directly, 
g = 2.26 JUt. 

As mentioned before, when the magnetic field exceeds 
the upper critical field, H C 2, the system is in the fully spin- 
polarized state (a corresponding dispersion of low-energy 
magnetic excitations is shown schematically in Fig. Q] (right)). 
Excitations from the ground state to the single-magnon state 
have been observed in our experiments and are denoted by 
squares in Fig. [2] (the resonance C). These excitations cor- 
respond to a single-spin flip from the S z = — 1 to S z = 
state and are uniformly delocalized over the entire lattice with 
a well-defined momentum k. The S z = Q state propagates 
along the lattice as a free quasiparticle with hopping J„ along 
the v direction (v = {a, b, c}), that arises from the transverse 
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FIG. 2: (color online) Top: Typical ESR transmission spectra in 
DTN, taken at frequencies of 216.9, 294.4, and 683.3 GHz at T = 
1.6 K, H || c. Bottom: Frequency-field dependence of magnetic ex- 
citations in DTN, taken at temperatures of 1.6 K (the modes A, B, C, 
and E) and 4.3 K (the mode F), H \\ c. Symbols denote the exper- 
imental results, and lines correspond to results of calculations (see 
text for details). 
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FIG. 3: (color online) Temperature-field diagram of the AFM- 
ordered phase obtained from magnetocaloric-effect measurements 
(opened squares), and the magnetization data taken at T — 16 mK 
(opened circles), H || c. Results of the quantum Monte Carlo calcu- 
lations of the phase-diagram boundaries and the magnetization using 
the set of parameters obtained as described in text are denoted by 
closed squares and circles, respectively. 



part of the Heisenberg interaction. There are also diagonal en- 
ergy gains of — 2( J c + 2J a ) due to the Ising part, and — D due 
to the single-ion anisotropy. The diagonal energy cost comes 
from the Zeeman interaction gfisH. Then, the single-magnon 
excitation dispersion can be calculated exactly: 



co{k) = gf i B H -D- 2(J C 



2J a ) + 2(J c cosk z 4 
J a cosk x + Jftcosfcj,) 



(2) 



The ESR transitions taking place at k = have the frequency 
u>c = — D. The best fit of the ESR data denoted 

by squares in Fig. [2] reveals D = 8.9 K for the anisotropy 
constant. From the exact expression for H C 2, given as 1 17] 



H c2 = —(D + A^J v ), 



9VB 



(3) 



and using H C 2 = 12.6 T 11161 11711 we obtain J v = 
J c + 2J a = 2.557 K. The zero-field dispersion of magnetic ex- 
citations calculated using neutron-scattering data 11711 yields 
Ja/Jc = 0.082. Thus, in addition to the anisotropy constant 
D = 8.9 K, all three exchange parameters, J c = 2.2 K and 
J a = Jb = 0.18 K, can be calculated quite precisely. 

The phase boundary (obtained from magnetocaloric-effect 
measurements) and the field dependence of the magnetization 
at T = 16 mK were computed for the obtained set of param- 
eters by means of a quantum Monte Carlo simulations for a 
finite lattice of L 3 sites, L = 16. Fig. [3] shows a very good 
agreement between the calculated (closed symbols) and ex- 
perimental (opened symbols) data. 

In addition to ordinary single-magnon states and two- 
magnon continuum, the theory |9] predicts the existence of 
two-magnon bound states (sometimes referred to as single- 
ion bound states |5J). The physical picture of the two-magnon 



bound-state excitations corresponds to a double-spin-flip tran- 
sition from S z = — 1 to S z = +1. The transverse part of the 
Heisenberg term of TL mixes this state with the one that has 
a pair of S z =0 states. Since the diagonal energy differ- 
ence between these two states, 2D, is much bigger than J c 
and J a (associated with hopping in the c and a directions), 
the distance between the two S z = sites remains finite, 
giving rise to a two-magnon bound state. The two-magnon 
bound states appear to be a specific feature of anisotropic 
spin-1 Heisenberg systems. It is worth mentioning that the 
two-magnon bound states were already predicted in 1970 by 
Silberglitt and Torrance ^ for Heisenberg ferromagnets with 
single-ion anisotropy. Later on, this subject attracted a great 
deal of attention due to its potential relevance to the intrinsic 
localized spin modes in anisotropic ferromagnets I22j] and an- 
tiferromagnets 12311 . It was suggested |9] that the two-magnon 
bound states should make a distinct contribution to the excita- 
tion spectrum of S = 1 large-D AFM chains above the upper 
critical field H C 2 and that their effect can be unambiguously 
identified by ESR measurements. A signature of two-magnon 
bound states was obtained by means of high-field ESR in 
the spin-1 chain compound Ni(C2HgN2)2Ni(CN)4 (known as 
NENC) 12411 . A broad absorption was detected in the high- 
field spin-polarized phase. Based on analysis of the tempera- 
ture dependence of the ESR intensity, this feature was inter- 
preted as transitions from the single-magnon to two-magnon 
bound states. However, a closer look at ESR spectra in NENC 
uncovered a more complicated structure than predicted (in- 
cluding, for instance, two non-equivalent Ni sites and, most 
likely, a finite in-plane anisotropy), which made the accurate 
quantitative comparison with the theory difficult. 

Here we report on the first observation of transitions from 
the ground state to two-magnon bound state in a spin- 1 AFM 
chain system with strong easy-plane anisotropy in the high- 
field FSP phase. The corresponding excitations are denoted by 
triangles in Fig. 12 bottom (the resonance E, Fig. [2] (top)). The 
frequency-field dependence of the ground-state two-magnon 
bound-state excitations can be calculated exactly using the 
set of parameters obtained as described above. Results of 
corresponding calculations are shown in Fig. [2] (bottom) by 
the line E. One more resonance absorption was observed at 
higher temperatures (which indicates transitions within ex- 
cited states). The corresponding data obtained at T = 4.3 K 
are denoted in Fig. [2] (bottom) by stars. This ESR mode corre- 
sponds to transitions from the single-magnon to two-magnon 
bound states (Fig. Q] (right)), which occur at k = 0. The 
frequency-field dependence of these transitions can be calcu- 
lated, using the expression uip = uje—^c (where loe and ojc 
are the excitation frequencies for the modes E and C, respec- 
tively), and is denoted in Fig.|2]by the line F. In both cases, ex- 
cellent agreement with experimental data was achieved, which 
shows that our model actually works fairly well up to T ~ J. 

Let us now discuss the low-field quantum-paramagnetic 
phase. As mentioned before, in this phase the system has 
the M s = ground state and a gapped excitation spectrum 
formed by AM S — ±1 transitions. The obtained set of spin- 
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Hamiltonian parameters can be used to calculate the lower 
critical field, employing the spin-wave theoretical approach 



H, 



1 



9(J>B 



In 2 -4s 2 nJ2 J v 



(4) 



The parameters fi = 10.3 and s 2 
the self-consistent equations: 



0.92 can be obtained from 



Ik, 

LO k ' 
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(H = 0). It is 



where 7k = 2 ^2 J„ cos k v and = toj 
important to mention, that the disagreement between the cal- 
culated (H c i = 2 T) and experimental (H c i = 2.1 T) values 
is due to limitations of the spin- wave theory applied to quasi- 
one-dimensional spin systems. Using the spin-wave theory 
1 1811 and the set of parameters obtained as described above, 
the frequency-field dependence of the ESR transitions can be 
calculated using the expression |QT 
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Eq. (IT) predicts a value of H c i = 2 T, that is in fact in rela- 
tively good agreement with the experimental value H c i = 2.1 
T. The results of the calculations are presented by the lines 
A and B in Fig. [2] (bottom) together with experimental data 
denoted by circles. The agreement in this case is reason- 
ably good, considering that the quasi-one-dimensional nature 
of the system causes pronounced quantum fluctuations of the 
AFM-order parameter in the low-field QPM phase. 

In summary, a systematic ESR study of the magnetic exci- 
tations in DTN, an S — 1 Heisenberg AFM chain material in 
the large-D regime, has been presented. Investigation of the 
high-field magnon excitations allowed us to obtain a reliable 
set of the spin-Hamiltonian parameters, which was employed 
to calculate the ESR spectrum in a broad range of magnetic 
fields and frequencies. These values agree very well with 
the ones obtained from fitting the AFM-phase boundary and 
low-temperature magnetization of DTN with results of quan- 
tum Monte Carlo simulations, including both critical fields. 
The parameters were used to calculate the frequency-field de- 
pendence of two-magnon bound-state excitations, predicted 
by theory and observed in DTN for the first time. Excellent 
agreement between the theory and experiment was obtained. 
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